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Abstract 
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1 Introduction 



The integrability of two-dimensional lattice models with periodic boundary condition is 
a consequence of the Yang-Baxter equation[l,2], 

R 12 (u - v)R 13 (u)R 23 (v) = R 23 (v)R 13 (u)R 12 (u - v) (1) 

where the R-matrix is the Boltzmann weight of the two-dimensional vertex model. As 
usual, R 12 (u), R 13 (u) and R 23 (u) act in C n <g> C n <g> C n with R 12 (u) = R(u) <g> l,R 23 (u) = 
1 ® R(u), etc. 

During the last years, much more attention has been paid on the investigation of inte- 
grable systems with nontrivial boundary conditions, which was initiated by Cherednik[3] 
and Sklyanin[4]. They have introduced a systematic approach to handle the boundary 
problem in which the reflection equations appear. In addition with the Yang-Baxter equa- 
tion, the reflection equations ensure the integrability of open models. Using this approach, 
Sklyanin [4], Destri and de Vega [5] solved the spin- | XXZ model with general boundary 
conditions by generalizing quantum inverse scattering method. Under a particular choice 
of boundary conditions, the Hamiltonian is U q [sl{2)\ invariant [6]. In [4], Sklyanin as- 
sumed that the R-matrix is P— and T— symmetric. Furthermore, the R-matrix satisfies 
unitarity and cross- unit arity properties. Because only few models satisfy these proper- 
ties, Mezincescu and Nepomechie [7] extended Sklyanin's formalism to the PT-invariant 
systems. Thus, all trigonometric R matrices listed by Bazhanov [8] and Jimbo [9] can be 
related to 1-dimensional quantum spin chains in this formulism. Using the unitarity and 
cross-unitarity properties of Belavin's Z n elliptic R-matrix, we have constructed the open 
boundary transfer matrix with one parameter [10,11]. 

On the other hand, the study of open boundary conditions in 2-dimentional field theory 
is related to the Sine-Gorden, Affine Toda and O(N) Sigma models [12,13,14,15]. Sklyanin 
generalized the hamiltonian to the case nonlinear partial differential equations with local 
boundary conditions [15]. The reflection matrix is consistent with the integrability of the 
systems. 

Recently, Foerster and Karowski have used the nested Bethe ansatz method to find 
the eigenvalues and eigenvectors of the supersymmetric t-J model with open boundary 
conditions and proved its spl q (2, 1) invariance [16]. Gonzalez-Ruiz also solved this problem 
with the general diagonal solutions of the reflection equation [17]. The investigated model 
is a graded 15-vertex model characterized by two bosons and one fermion. De Vega 
and Gonzalez-Ruiz have also generalized the nested Bethe ansatz to the case of SU q {n)- 
invariant chains[18]. 

The graded veretx model was first proposed by Perk and Schultz [19]. In this model all 
variables take m + n different values and the weights favor ferroelectric or antiferroelectric 
configurations. In the references [20] and [21], the Bethe Ansatz equations and the exact 
free energy and excited expectrum of this model with periodic boundary condition are 
found. The finit size correction shows the central charge of the model being m + n — 1 
(replacing n in ref [20] by m + n — 1 ). In fact, the supersymmetric t-J model is a special 
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Perk-Schultz model (m = 2,n = 1). Under an appropiate boundary condition, the model 
enjoys beautiful structure as quantum group symmetry. This motives us to consider the 
general graded vertex model with open boundary condition. 

In this paper, we use the nested Bethe ansatz method to find the eigenvalues and eigen- 
vectors of the transfer matrix for a graded vertex model with open boundary conditions. 
The transfer matrix with fixed boundary conditions is proved to be SU q (n\m) invariant. 
When m = 1, the model reduces into the q-deformed version of the generalized supersym- 
metric t-J model with n components. The hamiltonian contains a spin hopping term, the 
nearest neighbour spin-spin interaction and the contribution of boundary magnetic fields 
(see equation (24) ). Now, we outline the contents of this paper. In sect. 2 we introduce 
the SU q (n\m) vertex model. We find the matrices K ± which define boundary conditions 
and nontrivial boundary terms in the hamiltonian. The relation between the transfer 
matrix and the hamiltonian of the generalized supersymmetric t-J model is also discussed 
as an example. Sect. 3 covers to the diagonalization and the energy spectrum of the model 
with open boundary conditions in frame work of the nested Bethe Ansatz. In sect. 4 we 
show that the vertex model is a realization of the quantum supergroup SU q (n\m). A 
proof that the transfer matrix with open boundary conditions is SU q (n\m) invariant is 
given. In sect. 5 the summary of our main results is presented and some further problems 
are discussed. The appendix contains some detailed calculation. 

2 The vertex model and integrable open boundary 
conditions 

Our starting point is a graded vertex model which was introduced by Perk and Schultz 
[19]. The thermodynamics of the model with periodic boundary condition was studied in 
[20,21]. Some interesting application of this model in quantum field theory was considered 
by Babelon, de Vega and Viallet [22]. 

The model is defined by vertex weights R(u), whose non-zero elements are 



sin(r] + (-l) £a w), 
sin(u)(-l) eaeb ,a ^ b 
sin(ri)e lusi9n{a - b) ,a^ b. 

(2) 
where 

£ f 0, a = 1, • • • ,771 

I 1, a = m + 1, • • • , m + n, 

i] is an anisotropy parameter, a, b are indices running from 1 to m + n. For convenience, 
we denote 



K» = 
RtM = 
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a(u) = sin(u + rj) , b(u) = sin(u), , , 

w(u) = sin(r) — u) , c±(u) = sin(rj)e ±lu . 

This model is an m + n state vertex model characterized by m bosons and n fermions. 
If m — 2, n — 1, this model reduces to the one studied by Foerster, Karowski [16] and 
Gonzalez-Ruiz [17]. The A m _ x vertex model studied by de Vega and Gonzalez-Ruiz [18] is 
just special case n — 0. When m = n = 2, we can get a new electronic strong interaction 
model which is a generalization of the model proposed by Essler, Korepin and Schoutens 
[22]. 

The R-matrix defined by (2,3) is a trigonometric solution of the Yang-Baxter equation 
(1). The local transition matrix which is the operator representation of Yang-Baxter 
equation is the (m + n) x (m + n) matrix L(u) satisfying the following equation [19,20,21]: 



Ri 2 (u - v)L 1 (u)L 2 (v) = L 2 (v)L 1 (u)R 12 (u - v) (5) 

where L\(u) = L(u) ® 1, L 2 {u) = 1 <S>L(u). The standard row-to-row monodromy matrix 
for anJVxiV square lattice is defined by 



T(u) = L N (u) ■ ■ ■ Lx(«) 

= Ron(u) • • ■ Roi(u), (6) 

Throughout the paper, L{u) is assumed to be in the fundamental representation. T(u) 
also fulfills the Yang-Baxter equation 

Rn(u - v^iu^v) = T 2 (v)T 1 (u)R 12 (u - v). (7) 

The operator T is an (m + n) x (m + n) matrix of the operators acting in the quantun 
spaced- 

We can see that the R-matrix does not satisfy the Sklyanin's P- and T-symmetry, but 
fulfills PT invariance 

Pi2Ri2{u)P l2 = R t ^(u). (8) 
It also obeys the unitarity and cross-unitarity properties 

Ri 2 (u)R 2 i(—u) = sin(u + r))sin(r) — u) • id, (9) 

R%(u)M!R%(-u - d^M^ 1 = -sin(u)sin(u + drj) ■ id. (10) 
where d = m — n and M is a (m + n) x (m + n) matrix 
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M bc = 5 bc M b , 

( e^ 2 ^" 1 ), b<m 
M b = { (11) 
[ ed e i2 ^ 2m - b ). b>m 

This can be verified by straightforward calculation. One can also verify that M is a 
symmetry matrix of the R-matrix: 



[M ® M,R 12 (u)] = (12) 

Now, we can use Mezincescue and Nepomechie's generalized formalism to construct 
integrable systems with open boundary conditions. In our case, the reflection equations 
take the following form [7] : 

R 12 (u - v)K^{u)R 21 {u + v)Kz(v) = K^(v)R 12 (u + v)K^(u)R 21 (u - v) (13) 

R 12 {-u + v)K+{u) tl M^R 21 {-u -v- d^MxK+iyf 2 
= X 2 + (v)* 2 M 1j R 12 (-m - v - d^M^K+iuf'R^i-u + v) (14) 

Obviously, there is an isomorphism between K + (u) and K~{u). 

<f):K-(u)^K + (u)=K-(-u-^-yM (15) 

Therefore, given a solution K~(u) of equation (13), we can also find a solution K + {u) 
of equation (14). But in a transfer matrix of an integrable lattice, K~{u) and K + (u) need 
not satisfy equation (15). In this paper, we will take equation (15) to define K + . 

After a long calculation, we find a solution of the reflection equation (13) 

K-{u) = id (16) 

Correspondingly, 

K + (u) = M (17) 
Taking Sklyanin's formalism, the double-row monodromy matrix is defined as: 

U{u) =T(u)K-(u)T- 1 (-u) (18) 

where T _1 (w) is the inverse of T(u) in the auxiliary and quantum spaces, which explicitly 
is: 
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T-\u) = L- l \u)---L- N \u) 

= R 01 1 (u)---R^(u) (19) 

With the help of the Yang-Baxter equation (7) and the reflection equation (13), one can 
prove that the double-row monodromy matrix satisfes the reflection equation 

Rn(u - v)U 1 (u)R 21 (u + v)U 2 (v) = U 2 (v)R 12 {u + v)U 1 (u)R 2 i(u - v) (20) 
In this case, the transfer matrix is defined as: 



t(v) = trK + {y)U{y). (21) 

Using the reflection equations (14,20) and the properties of the R-matrix (8-12), one can 
prove 

[t(u),t(v)} = 0. (22) 

So the transfer matrix constitutes a one-parameter commutative family which ensures 
the integrability of the model. As indicated by Sklyanin, the transfer matrix is related to 
the hamiltonian of the quantum chain with nearest neighbour interaction and boundary 
terms 



1/(0) = 2trK + (o)—^—H - 2Ncotr] (23) 
sin(rj) 

From equation (21), one can derive the explicit expression of the hamiltonian, which is 
omitted here because it is not used in the following discussion. In order to compare it with 
the SU q (2\l) supersymmetric t-J model, we give the hamiltonian under m — 1, which is 
defined: 



N-l 

j=l s 

N-l / rn 

+ J2 COS7 l Yl n han j+ l,a + COST] ■ (rij + U j+1 ) - COST] ■ (UjUj+i) 
j=l \a=l 

N-l I 

3=1 \«€A + 
N-l 

+ E isinr l ■ (E n ho,n j+ i, h - E n jta n j+ltb ), (24) 

j=l a<b a>b 
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where cj s (cj S ) creates (annihilates) an electron with spin component s,s — 1, 2, • • • , n+m— 
1 located j-th site, rij is the density operator, Sj is the spin matrix at site j. A + denotes 
the set of positive roots of the su(m) algebra. P is a operator projecting out doubly 
occupied states. The constraint is that more than one electron on each site is strictly 
prohibited. As we know [16,17], this hamiltonian is not hermitean, but it possesses real 
eigenvalues. We will show the hamiltonian to be SU q (n\m) invariant. 



3 Nested Bethe ansatz for open boundary conditions 

The graded vertex model with periodic boundary condition was investigated by de Vega 
and Lopes [20,21]. Based upon the Yang-Baxter equation, they obtain the Bethe Ansatz 
equations by using the nested Bethe ansatz method (periodic case). In this section, we 
want to generalize the nested Bethe ansatz method to solve the eigenvalue problem of the 
transfer matrix (21). In this case, the operator commutative relations are ruled by the 
reflection equation instead of the Yang-Baxter equation. 

As we know, the double-row monodromy matrix satisfies the reflection equation. It 
is convenient to denote tt_ = u — v , u + = u + v . We rewrite the equation (20) in the 
component form: 

i?i2(«-)^ 2 2 f/(«) cl d lJ R2i(« + )S^(t;), 2fe2 

= f/(^)a 2C2J Rl2(« + )^ 2 2 f/(«) cl(il i? 2 l(«-)K 2 2 (25) 

where the repeated indices sum over 1 to m + n. Next, we introduce a set of notations 
for convenience: 



A(v) 


= U(v] 


In, 


B a (v) 


= U{v) 


•la, 


C a (v) 


= U(v] 


lal, 


D ab (v) 


= U(v] 


lab, 



From equation (25) we will find the commutation relations. In order to simplify these 
relations, we introduce new operators: 



D ab (v) = D ab (v) - S ab ff V) ^ A(v) (27) 

Considering the vertex model defined by equations (2,3), we rewrite equation (27) in an 
explicit form: 
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D ab (v) 



D ab (v)-5 ab ^j^A(v), m^O. 
a(2v) 

D ab (v)-5 ab ^lA(v), m = 0. 
w(2v) 



(28) 



After some tedious calculation, we have found the commutation relations between A(v), 
D ab (v) and B a (u) (a, b — 2, • • • , m + n). The final results take the form (see Appendix A) 



A(v)B b (u) 



a(u — v)b(u + v) 
a(u + v)b(u — v) 

c+{u — v)b{2u) 
a(2u)b(u — v) 

c_(w + v) 



B b {u)A{v) 
B b (v)A(u) 



a(u + v) 



B c (v)Dd,(u), 



(29) 



+ M2,. 1 + ,)^!g? ^' 2 " + ^B^m. (30) 

All indices take values fron 2 to m + n, and the repeated indices sum over 2 to m + n. 
The commutation relations presented above are only applicable to the cases m > 1. If 
m=0, the commutation relations change to the following form: 

w n ^ / \ w(u — v)b(u + v ) „ . , , . , 
u>(w + v)b(u — v) 

c + (u — v)b(2u) „ . . . . . 

y -5 b (t;)A(M) 



w(2u)b(u — v) 
C-(u + v) 



w(u + v) 



B c (v)D cb (u) (31) 
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The rule for indices is the same as the one in equations (29,30). 

It is easy to find the so-called local vacuum ef. We call the direct product of local 
vacuum a reference state or vacuum state. It takes the form: 



|m C >=n(l,0,--,0) ( , (33) 
where t denotes the transposition. One can find 



A{u)\vac> = a(u)\vac>, 

C a {u)\vac > = 0, 

B a {u)\vac > 7^ 0, 

a(u) = [R(u)[\] N [R-\-u)[\] N . (34) 

Next, let us calculate the action of D ab (u) on the vacuum state. We first recall the 
definition of D ab (u), and find 

D ab (u)\vac >= T(u) al T-\-u) lb \vac > +T(u) ac T~\-u) cb \vac > . (35) 

The contribution of the first term can not be calculated directly. We will use the following 
method to find it. Taking v = —u in the Yang-Baxter equation, we can get: 

T 2 - 1 (- M )i? 12 (2 M )T 1 ( M ) = T 1 ( M ) J R 12 (2 M )T 2 ~ 1 (- M ) (36) 

Taking special indices in this relation and applying both sides of this relation to the 
vacuum state, we find: 

T(M) a iT _1 (-M)i fe |wac >= ~t~t~tt (8 ab a(u) - T(u) ac T _1 (-u) c6 ) \vac > . (37) 

K(2u) u v / 

Substituting this relation to eq.(35), we have the result: 

D ab (u)\vac >= M|^r«W + (l - l^yrr) b N (u)b N (-u)}\vac > . (38) 
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So we have 



D ab (u)\vac >= 5 ab - ^y|) b N {u)b N {-u)\vac >, (39) 

where b(u) = R~ l (u) a ^ h , a ^ b. In conclusion, the results of the action of A, B a , C a and D ab 
on the vacuum state are listed as: 

A(u)\vac> = [R(u) 1 1 1 1 ] N [R- 1 {-u)\ 1 1 } N \vac>= a{u)\vac> 
( c + (2u) \ 

D ab (u)\vac> = 5 ab f 1 - j^ u yA b N (u)b N (-u)\vac >= 5 ab (3(u)\vac > 
C a (u)\vac > = 

B a (u)\vac> ^ 0. (40) 

Note that the action of B a (u) on the vacuum state is not proportional to the vacuum 
state. 

We show that the eigenvectors of transfer matrix t(u) can be constructed by repeatedly 
applying operators B b .(vi) on the vacuum state 

• • • , v L ) = B bl ( Vl ) ■■■B hL (v L ) \vac > F bl "' bL (41) 

Before using the Bethe ansatz method, let us introduce a set of notatons that will be used 
in the following. We denote 

R2i(u)? d /R(u){ 1 1 = R 21 (u)? d . (42) 
So from Appendix A, the commutation relations between B's take the form 

B bl ( Ul )B b2 (u 2 ) = R 12 ( Ul - u 2 ) d b ^B d2 (u 2 )B dl ( Ul ). (43) 

By repeatedly using this relation, we can commute B{vk) with B(vk-i), • • • , B(vi), re- 
spectively. 

B bl (vi) ■ ■ ■ B bL (v L )\vac > 
= S(v k , {ui}) J B dl {v k )B d2 (ux) • ..B dk (v k . 1 )B dk+1 (v k+1 ) ■ ■ ■ B dL (v L )\vac > (44) 

where 



S(v k , {v t })tX 

= II ^ jdj RMvi - v k ) d c £R 12 (v 2 - v k ) c c f 2 ■ ■ ■ R 12 (v k ^ - v k ) c b k kb l d l (45) 

j=l+k 
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Here d\ and b k are considered as the "auxiliary space" indices, &i, • • • , fcfc-i, fcfe+i, • • • , &l 
and c?2, • • • , di are the "quantum space" indices. Notice that in the six vertex model, the 
B's are commutable with each other. So one can use the symmetric argument that v k and 
Vi are equivalent to each other. Now, we know that the B's are not commutable with each 
other, but the relation (44) ensures that we can also use something like the symmetric 
argument. Actually, we can see from this relation that v k and v\ are in an equivalent 
position if we omit the function S. 

In the following, we deal with the case of R(u)\\ = a(u). It is convenient to introduce 
the notation: 



~ ,, \ a h, = R2l(-Vl - v ijg 2 b 2 

a2b2 a(vi + v i )a(-v 1 -ViY 



R12(V1 + V L ) C q L L - p l dL R2l(vi ~ V L )\Z\bAM) 

— [^(fii.^^^o^H-fii,^)- 1 ];^. (46) 



a(fi-Wi-vi) 9L_1 L 

Here we have used the unitarity properties of R matrix, and %i = v j + r]/2, {fl{v\)) ah = 
S a bP(vi). 

Now, let us evaluate the action of A(u) on Following the algebraic Bethe ansatz 
method, many terms will appear when we move A(u) from the left hand side to the right 
hand side of B a 's. They can be classified in two types: wanted and unwanted terms. The 
wanted terms in A(u)^/ can be obtained by repeatedly using the first term in relation 
(29), the unwanted terms arise from the second and third terms in relation (29), they are 
the types that v k is replaced by u. One unwanted term where B(v\) is replaced by B{u) 
can be obtained by using first the second and third terms in relation (29), then repeatedly 
using the first terms in relation (29) and (30). Using this results we can obtain the general 
unwanted term where B(vk) is replaced by B(u). So we can find the action of A(u) on ^ 



A(u)B bl ( Vl ) ■ ■■B bL (v L )\vac > F bl '" bh 
-A- a(v ~ — u)b(v ~ + u) . . „ . . „ . . h 

= jg ^; -«.)«!»; +U ) a(u) • BM • • ■ BMae > 

! j^ -c+jvk - u)b(2v k ) -p. ajvj - v^bjvj + v k) a ^ F \n-b L 
fc=1 a(2v k )b(v k - u) j=1 ^ k a(vj + v k )b(vj - v k ) k 
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•B a (u)B a (bi) ■ • ■B dl (v l -,)B it+1 {v t+1 ) ■ ■ ■ B dL (v L )\vac > S(v k , {v,})^ 



n ^|pL fe)%,w);:;: (47) 

•^(-^l) -1 • • .L^\-%rv k -i)- 1 LV{-v k ,v k+ i)- 1 • • • f) L )- 1 ]^ ;£ 

(«i) • • •S djfc (u fc _ 1 )S dt+1 ( Ufc+1 ) • • .fl dj > L )|vac > F bl -^ 
Recalling the definition of the transfer matrix, we rewrite the transfer matrix as: 

m+n 

t( U ) = Y, K a(u)U(u) aa 

£1=1 

m+n _ (m+n Icy \ ^ 

= £ Kt{u)D aa {u) + + £ K aH^4 + 

a=2 I a=2 \ / ) 

= J? KMDUu) + sin ^ + d ^ e ^A(u) (48) 
^ sin{2u + 77) 

Before calculating the action of the transfer matrix on we should evaluate the action 
of K+(u)D aa (u) on it, which reads 

m+n 

K+(u)D aa ( U )B bl ( Vl ) ■ ■ ■ B bL (v L )\vac > F b ^ 

a=2 

m+n L 1 

= V K+(u)F h - bL f3(u) TT 77 T77 r 

h , f} 1 b(u-v j )b(u + v ] + V ) 

■Ri 2 (u + v 1 + V ) a P t R ^ u - «i)STi2i2(« + v 2 + v) P P l£R2i(u - v 2 ) s X 

■ ■■r 12 (u+v l + v y p ^R 21 ( u - v L y s i% L 

B dl (v i)B d2 (v 2 ) •■•B dL (v L )5 PnSn \vac > +u.t. (49) 

where u.t. stands for unwanted term. Using the definition of lP-\u,Vj) and its inverse, 
we rewrite relation (49) as: 

h aK fAKu-Vj^u + Vj + rj)^ 



bi-b L 

B dl (v^B^ (v 2 ) ---B dL (v L ) \vac > +u.t. (50) 



Here 



13 



di—d L 
b!-b L 



u, tfx) - - - (fi, £ L )£ (1) (-«, ^l)- 1 • • • L^(-u, v 1 )- 1 Y 1 ''T } ( 51 ) 

As mentioned above a is the index of the auxiliary space and bi, are the indices of the 
quantum space. The unwanted terms in equation (49) take two forms. After some long 
tedious calculation based upon the similar considerations as in the A(u) case, we can get 
the following expression After long tedious calculation 

u.t. = ELi^K,{^})tt iF6l ''' 6i -«K) e ^ 

c + (u + v k )b(2v k )b(2u + drj) a(vj - Vk)b{vj + v k ) 

a(u + v k )a(2v k )b(2u + rj) a(vj + v k )b(vj - v k ) 

■B dl (u)B d2 (vi) ■ ■ ■ B dk (v k -i)B dk+1 (v k+1 ) ■ ■ ■ B dL (v L )\vac > 
+ ELi S ( u k, { v i})tl'.. C b L L Fbl ''' bL ■ P(v k )e id ^ 

-c + (u - v k )b(2u + dr]) a(vj - v k )a(v k - Vj) (52) 

b(u - v k )b(2u + rj) b(v k - Vj)b{vj + v k + rj) 

(Vh) • ■■LW(v k ,v k - 1 )Ll 1 ^v k ,v k+1 ) ■ ■ ■ 
■ ■ ■ LW(v k , v L )L^(-v k , v L )~ l ■ ■ ■ v^)- 1 

Lm(-v k ,v k _ 1 )- 1 ...Lm(-v k ,v 1 )- 1 " 



B dl (u)B d2 (yi) -"B dk (y k -i)B dk+1 (v k+1 ) ■ ■ ■ B dL (v L ) \vac > . 
In order to simplify equations (47) and (52), we need an important relation 

{L« (v k , (Vh) • • • L« (v k , vh-JlW (v k , v k+1 ) L« (v k , v L ) 

■L^i-v^VL)- 1 ■ ■■L^(-v k ,v k+1 )- 1 L^ 1 \-v k ,v k _ 1 )- 1 ■ ■ ■ 

> Ci—Cl 

• S/,,(,/) :e-^5^,{^})-;-r( 2 )(^,{^})-;^ (53) 



where 



sin(2v k + dr)f ^^>h-H' WW)* 



r {2 \u,{v l })ti:t L 

m+n 

= Y, K a(y){{L {l \uM)---L^\u,v L ) 

a=2 
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.L«(-«,iJ L )- 1 ...L( 1 )(-«,tJ 1 )-X7} 

E^(«){(r (1 H«,W)r«(-«,{i; i })-^7} (54) 



Using the equations (48) and (53), we then obtain the action of t(-u) on ^ 



v ; v ; sin(2u + r]) jJib{vj -u)a{vj + u) 

F b ^B bl ( Vl )---B bL (v L )\vac> 

+ n / vj "w o( r v i\ ^ r(2) ^ 

/=i o( u — v j)b\ v j + u + rj) 
F bl "^B dl (v 1 )---B dL {v L )\vac> 

+ y f -c+(vk ~u) c _ iv + c+(v k + u) \ sin(2u + dri)b(2v k ) cidn 
fr{ V K^fc - u ) a ( v k + u) j sin(2u + r])a(2v k ) 

■ n p'-^'l^U vtWvtoMKXf** (55) 

i= i^ fe 6(uj - v k )a(vj + v k ) 

■B dl (u)b d2 (ui) • • • 5 dfe (w fc _i) J B dfc+1 (ufc+i) • • • -B di (u L ) |vac > 
_ A / c_(m + ^) _ i7? c+(u- v k )\ sin(2u + fir?) 



fc=i 



a(w + i>fc) 6(m — t>fc) y sin(2u + rj) 



sin(rj) j^- a(v k - Vj)a(vj - 



sin(2v k + rj) b(v k - + v fc + v) 

B dl (u)b d2 (vi) ■ • • B dk (v k -i)B dk+1 (v k+ i) ■ ■■ B dL (v L )\ vac > 

From the above equation, one can see that the function \I/ is not the eigenstate of t(u) 
unless F's are the eigenstates of and the sum of the third and the fourth term in the 
above equation is zero, which will give a restriction on the L spectrum parameters 
So, we have the following results: 

If F is the eigenstate of with the eigenvalue satisfmg equation (57), then \1> 
is the eigenstate of t(u) with the eigenvalue A^, 

A d) M = a ( u Y(d~Dv sl < 2u + dr i) A a fa - u ) b ( v J + u ) 

sin(2u + rj) b{ Vi - u)a{ Vj + u) 



,=1 K M _ V j)KVj +U + T})' 
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where 

A^(u, {v t })F 

a(v k )b(2v k )sin(2v k + drj) i(d _ 1)v bjvj + v k ) (57) 

f3(v k )a(2v k )sin(r)) a(v k - vj) 

Therefore, the diagonalization of t(u) is reduced to finding the eigenvalue of t^ 2 \ The 
explicit expression of (see equation (55)) implies that can be considered as the 
transfer matrix of an L-sites quantum chain, in which every spin takes m + n — 1 values. 
The related Yang-Baxter equation is the same as the one of t(u), exception R being an 
(m + n — l) 2 x (m + n — l) 2 matrix. Hence, we can use the same method to find the 
eigenvalue of t^ 2 \ Repeating the procedure m times, one can reduce to a subsystem r^ m+1 ) 
which is an n x n matrix in auxiliary space. The related Yang-Baxter equation is also 
defined by equation (2), but one should notice that in this case all e a = — 1 due torn = 0. 
In order to diagonalize 

r (m+i) j we 

need the definition of D by the second equation (28). 
The elements of satisfy equations (31) and (32). Following the same procedure, one 
can further reduce the r^ m+1 ^ into the r*™ 4 " 2 ) subsystem. The late has the same structure 
as the former. In this case one finaly obtains the eigenvalue of T ( m + n - 1 ) . This is the well- 
known nested Bethe Ansatz. Because the wave-functions are not needed in this paper, 
we omit them here. The eigenvalue and the constraint on the spectral parameters read as 



r«M^})F = 

a (2) K-,R» = 



AW(u,{vt 1] },{vl k) }) 



V ' X 1 i! sin(2u + kr]) 

g a{vf ] - u)b{u + vf ] + (k - 1)77) 



L i b(vf ] - u)a(u + v) K> + (k — l)r]) 

5> j=\ b(u - v\ k) )b(u + vf ] + krj) 
■M k+1 \u,{v\ k) },{vl k+1) }) 

(1 < k < m) 



(k) 



(58) 



16 



sm(—2u + (k — 2m)r]) 

g ) - «)6(u + ) + (2m - fc + 1)77) 
j=i b(vf^ — u)w(u + Vj k ^ + (2m — k + 1)77) 



=i b(u — Vj )b(u + Vj + (2m — fc)?7) 
■A ik+1 \u,{4%{vt 1] }) 

(m + 1 < k < n) 

(59) 



and 



A (*+!)/.,(*) r.,(*h r 7J (*+i)iN . « (fc) (4 fc) ; fa" sm(2?j fc) + dq) i{d+k _ 2)r) 
{Vl ' {Vi 1{V * !) " ^(^fa^}) 

swgtf + (fc - 1)77) ^ sm^gj + ?/f } + (fc - 1)77) 

sin(2vi^ + krf) j=i,^i sin(vj k ^ — — 77) 
(1 < fc < m) 

(60) 



A (/c+l)/ ? .(fc) r (fc)-. r ? ,(*+l)lN _ ^ 'fa » i(d+2m-fc)q 

>fa I) 

sin{2v\ k) + dri)sin{2v[ k) + (2m - fc + 1)77) 

sin{ri)sin{2v\ k) + (2m - fc)?;) ( 61 ) 
^| sin(u, (fc } + } + (2m — k + 1)77) 

11 • / (fc) (fc) . X 

j=i,# sin[v\ - v) + r/) 

(m+l<fc<m + n — 1) 
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where 



simu 



.(*-!) 



+ £77) 



II 

j=i sinyUj — u + rj) 



' ' 1 sin(w + vf ^ + (2m — k)rj) 



n 



sin(vf ^ — u — rj) 



1 < k < m 



m + 1 < k < n 



p t T 1 siniu + vf 1 - ) + (k — l)rj) sin(u — vf 1 h 
II 



sin[v 



(fc-i) 



u + rj)sin(u 



.(*-!) 



' + 77) 



sin(2u + (k — l))rj) 



sin(2u + krj) 



-IT] 



Pfc-1 

n 



sm(u + fj fc ^ + (2m — fc + \)rj)sin(u — vf 1> ) 



1 < k < m 



sin[Vj 



u + rj)sin(u — vf ^+77) 



sira(2u + (2m - fc + l))r?) 
sin(2u + (2m — fc)^) 



-177 



m + 1 < fc < n 



(62) 



In the above representation, we have assumed fj 1 ^ = Vj,vf^ = 0, P Q = N,Pi 



L and 



^(m+n+i) _ ^ Notice that (3^ k \u, {vf 1 }) vanishes at the special points v\ k ^ due to the 
factor sin(u — v\ ~ ') appearing in fr k \ Taking u = vf ^ in formulae (58) and (59), we 
can get another kind of constraints on 



.(*) 



(fe-1) rffc-l)! 



t sin(2vf ^ + <ir/) 
sin(2v\ k ^ + ^77) 



p t • / (fc) 



+ rj)sin(vf^ + i>}" ±J + (A; — 1)77) 



(fc-i) 



n 

1 < k < m 



(k) (k-l)\ ■ I (k) . (fe-1) . , \ 

f z )stn{Vj + v I + krj) 



simv 



(63) 



AW^-D, {w (*-i) }> = -aW^-D, {v <*-i) })e i<*w»-*)., + ^ 

sin(2v\ + (2m — fc)?7) 

j^. sin(Vj k ^ — v\ k ~ 1 ^ — r])sin(vf^ + v\ k ~ 1 ^ + (2m — k + 1)77) 

j=i sin(v^p — 1 ')sin(?) J -' ! ' + ^ + (2m — fc)?7) 

m + l</c<m + n 

(64) 

Now, changing the index k into + 1 in the above formulae, we can obtain constrains 
on A^ +1 \ Comparing these with equations (60) and (61), one can derive out the following 
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Bethe ansatz equations 



sin{v\ k ^ — if ^ + i])sin(vi k ^ + if ^ + krj) 



(fc) 



(k-iy 



II . , fW a-iK . , (fc) + Vj k ~^ + (fc - 1)77) 
"^i 1 sm(if -* — Vj k+lS) — rj)sin(v\ k ^ + i>f + krf) 
' M sm{v\ k) - vf +l) )sm{v\ k) + +1) + (k + 1)77) 
g sin(v\ k ^ — if ^ — r])sin(v\ k ^ + if ^ + (fc — 1)77) 
j=i,# sin(v\ k ^ - vf^ + r])sin(v\ k) + i;f ^ + (fc + 1)77) 
(1 < fc < m - 1) 

•Ji™ 1 ^ — i> ? - m ^ + i])sin(vi m ^ + i> ? - m 1 ' ) + mi]) 



n 



simvi 



sin{v\ 



(m) 



(m— 1)\ • / (m) (m— 1) 

7jj )sin{vl +Vj 



(m - 1)77) 



+ 1 sm(i>f^ — 7j- m+1 ^)sm(7j^ + i>- m+1 ' ) + (m — 1)77) 



J=i sin(v\ m ^ — t>f 1+1 - ) + rj)sin{v\ m ^ + vf 14 " 1 - 1 + 777,77) 



(65) 



(66) 



sin[v 



r k -\ 

n 



(if ^ — if ^ — r])sin(vi k ^ + i>f 1 ' ) + (2m — fc)?7) 



(fc) 



sin[v 



(k) 
1 



+ v 



(fc-i) 



sin[v 



n 

= n 



(fc) 
/ 



+ (2m - fc + 1)77) 
if +1 ^)sin(v\ k ^ + if + (2m — k — 1)77) 



(fc) 



(fc+i) ^ ^ s i n [y^ 4- t)^ fc+1 ^ + (2m — fc)?7) 



sin [v 



(fc) 
1 



.(*) 



i])sin(v^ + v { - K> + (2m — fc — 1)77) 



,(*) 



• / (fc) 



.(*) 



+ rj)sin(v\ k) + uf } + (2m - fc + 1)77) 



(67) 



(m+l<fc<m + n— 1) 

The above Bethe ansatz equations are very complicted, but they can be simplified by 
introducing the following new variables 



(fc) 
v) = < 



w 



(fc) 



£77/2, 



1 < fc < m 



w 



(fc) 



(68) 



(2m — k)rj/2, m + l<k<m + n 



The Bethe Ansatz equations then take the form 
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P L [ 1 sin(v\ kS> — vf 1 - > — r]/2)sin{v\ kS> + vf ^ — rj/2) 
/= i sin{v (k) - vf~ l) + /2rj)sin{v\ k) + vf~ l) + rj/2) 
-i sin ( v ( fc ) _ V f+V - rj/2)sin(vl k) + vf +l) - rj/2) 



n 



i sm(u, (fc) - vf +1) + r)/2)sin{v\ k) + +1) + 77/2) (69) 



Pt • / (*0 (fe) \ • / (fe) , (fe) \ 

sinyvi — Vj — rj)svn[vl + v) — rj) 



■ ( (k) (k) . X • / (fe) . (fc) . X 

j =1 ^1 sin\yl v) + r))sin(Vi +v^ + r]) 
(1 < k < m + n — 1, k ^ m) 

p - 1 sm(v\ m) - v\ m - 1] - V /2)stn(vl m) + vf^_ - V /2) 
j=i sin(v\ m ^ — f J - m_1 ' ) + /2r])sin(vl m ^ + t>j m ~ 1 - ) +r]/2) 

' Ik sin{v\ m) - vl m+1) - r]/2)sin{v\ m) + v^ m+1) - 77/2) 
= 1 

The function AW(u, • • •) must not be singular at u — vf^ (1 < j < p k , 1 < k < m + n — 1) 
since the transfer matrix t(u) is an analytic function of u. In fact, the equation (57) comes 
from the condition under which the unwanted term vanishes. One can understand this 
constraint from another point of view: From equation (56), we know that u = Vj = v^p is 
a pole of A^\u). In order to keep the analyticity of A^(u), one should need the residuce 
of A^(u) at Vj vanishing, which also gives the constraint (57). So, A^(u) is analytic 
at Vj. Similarly, equations (69) ensure the analyticity of A^\u) at all vf\ Therefore, 
the eigenvalues of the transfer matrix are analytic functions if the previous Bethe ansatz 
equations are satisfied. 

The energy spectrum of 1-dimensional quantum system defined by equation (23) can 
be derived fron A^(u). It is 

E - P y 1 1 M(d-m 

^ cos(rj + 2v k ) — cos(rj) Asin 3 (rj) 

In order to comparing our results with the Bethe Ansatz equations given in references 
[20,21], we introduce new varibles 

l<J<P k 
Pk + l<J< 2P k 

Then equations (69) and (70) deduce to the Bethe Ansatz equations (see, for example, 
equation (4) in ref. [20]) up to a phase. In this sense, the Bethe Ansatz equations for 
the system with quantum group symmetry is the duble of the ones for the same system 
with periodic boundary condition. One should note the constrain in the right hand side 
of equation (69), which will contribute a non-zero term to the free energy. 
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4 Quantum group structure of the model 



In this section we will show that the vertex model under consideration is a realization of 
quantum supergroup SU q (n\m), and we will also prove that the transfer matrix for open 
boundary conditions is SU q (n\m) invariant. 

Firstly, denoting x = e tv , q = e tr >, the Yang-Baxter equation becomes 

Ruix/yfr^TM = T 2 {y)Ti{x)R 12 {x/y). (72) 
We write the R-matrix as 

R(x) = xR+ - x~ x R-, (73) 
similaryly, the L operators can be written as 

L(x) =xL + -x~ x L-. (74) 
From the definition of R-matrix, L± can be written in the following form. 



L 1 = I qW% _ 1 ~ m ' (75) 
+l 1 (Tiq Wi , m < i < m + n, 



(?- 




-l) ? -5Ej 


(?- 


q~ 


-l) g -2Ej 


(?- 


q~ 


_1 )g2 Ej^. 



i < m, 
i = m, 

m < i < m + n, 



(76) 



i < m, 
i > m, 



(77) 



i—(q — q x )&iq 2 Wj , i < m, 

-(g-g~ 1 )e m g^ E ^-+ lWj+ "' m+1 , ^ = m (78) 
— (g — g~ 1 )crjejg~ 5 ^-'^M+ 1 ™', i > m. 

Here L ± are lower and upper triangular matrices with L + l - = LJ { — 0, if i < j, w iy i = 
1, • • • , m + n; e^, fi, i = 1, • • • , m + n — 1, are the generators of the SU(n\m) superalgebra 
in the graded Cartan-Chevalley basis; the definition of the matrices <7j and the details of 
the classical simple Lie algebra SU(n\m) are given in Appendix B. Recall the definition 
of the monodromy matrix T(x): In the limit x — > oo, 0, we find the leading terms T± of 
the monodromy matrix T{x) to take the form, 
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T. 



q -N/2 gW% 



i < m, 



Oiq Wi , m < i < m + n, 



(79) 



rp i+1 



a^q 2^j^i,i+i Wj Fi, % < to, 



a-q 

{ q- N a + q*^i* i > i + lVl 1 ~ 



i > m, 



(80) 



^J 1 * — 



q N / 2 (Jiq Wi , m < i < m + n, 



(81) 



T_i +1 = <| -^^i^ 1 , i = m, (82) 
—q N a-cr i E i q^ 1 ^^•.<+ 1 Wj , m < i < m + n. 

Here T± are lower and upper triangular matrices with T + * = = 0, if % < j, a± — 
q ±l l 2 (q — q-i), and 



o; t = (Ti (g fTj (8) • • • Cg <jj, i — to + 1, • • • , to + n, 

q ±W* = g ± Wl ^_^ g ±w^ i = l } ... >m _|_ n} 

AT 

Xj = ^ g^ 2 <g • • • <g g^ l/2 <g rrf ^ ® q h ' /2 ® ■ ■ ■ ® q h ' /2 ,i < to, (83) 
j'=i 

AT 

X m = Y,^ hm/2 ®---®^~ hm/2 ® X t®( a m+ iq hm/2 )®---^^rn+iq hm/2 ), 
J'=l 

AT 

*i = E ^ l/2 ® • • • ® ^' /2 ® ® (<Wi<W2 • • • <7 i+ ig _hi/2 ) 

J'=l 

(g • • • (g (o- m+ io- m+2 • • ■ o"i+ig _/ll/2 ), m <i <m + n. 

where X; = E u F h x { = e h , respectively. hi = Wi- w i+1 , % ^ m, and h m = w m + w m+1 . 
In the case of N — 2, these formulae define the coproduct of a Hopf algebra. From this 
point of view, the equation (83) can be written as 



X t = A (xi) = (A <g id)A N ~ 2 (xi) (84) 

In the following we will discuss the algebraic relations of (q Wi ,Xi). Taking the appro- 
priate limits of the R-matrix and the row-to-row monodromy matrix T, we have 
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lim xR(x/y) = —yR~ (85) 
lim x^Rix/y) = -R+ (86) 



x^oo 



y 



and 



\imx N T(x) = -T_ (87) 
Jim x~ N T(x) = T+ (88) 



In the limits x — > 0, oo, the Yang-Baxter equation gives: 

R±T 1± T 2 (y) = T 2 (y)T l± R ± (89) 

and 

i? ± T 1± T 2e = T 2£ T 1± R ± (90) 

with £ = {+,—}• These spectral-parameter-indepedent Yang-Baxter relations govern 
q-(anti) commutation rules and q-Serre relations for the quantum supergroup SU q {n\m). 
Substituting the definition of R±,T± into equation (90), we get 

q H iq H i = q H iq H \ 
q Hi Fjq- Hi = q a ^Fj, 
q Hl E jq - Hl = q~ aiJ Ej, 

q Hi q H i 

[Fi, Ei] = — , % ^ to, 

q-q 1 



q-q' 1 



E 2 m = F 2 m = ^[F l ,E 3 ] = 0,i^j, 
(FtfF^-iq + q-^FF^F + F^iFtf = 0, 
(EifEi^-iq-q-^EiEi^Ei + Ei^Ei) 2 = 0. 



(91) 



where Hi = Wi — W i+1 , Hi ^ to, H rn = W m + W m+ i, and is a component of the Cartan 
matrix which is given in Appendix B. The generators H^E^Fi, i = 1, • • • , m + n — 1, 
and relations listed above provide a definition of the quantum supergroup SU q (n\m). In 
the remaining part of this section, we will verify that the transfer matrix t(y) with open 
boundary conditions is SU q (n\m) invariant. The entries of lower and uper triangular 
matrix T± are elements of SU q {n\m). So, it is not necessary to compute commutators of 
t(y) with individual SU q (n\m) generators. If the relation 

[t(y),T ± ] = (92) 



23 



is correct, we are led to the conclusion that the transfer matrix t(u) is SU q (n\m) invariant. 
From eq. (89) we have the result 



[R ± T 1± ,T 2 (y)T 2 -\y- 1 )]=0 (93) 
Recall the relation (12), we have 

[R±, MiM 2 ] = (94) 

Similary, from the unitarity and cross-unitarity relations (9,10), with the help of PT 
invariance of R-matrix, we find 



R ± R^ = 1 
R^MiR^Mi 1 = 1 (95) 

So we have the identity = (-R^ 1 )* 2 which implies that the following relation 

M^^ 1 )' 2 ^^ = 1 (96) 

is correct. Notice that we choose K_ = 1 in this paper, so the transfer matrix can be 
written as t(y) = trMT(y)T^ l (y^ 1 ). Now, let us prove relation (92) 



T 1± t(y) = tr 2 T 1± M 2 T 2 (y)T- 1 (y- 1 ) 

= tr 2 M 2 R- ± l R ± T l± T 2 {y)T^\y- 1 ) (97) 

here we have added an identity R~±R± in the relation, then using the relations (95) and 
(96), we find 



••• = tr 2 Mi 1 RZ 1 M 1 (M 2 T 2 (y)T- 1 (y- 1 ))R±T 1± 

= tr 2 {Mr 1 J Ri 1 M 1 }' 2 {(M 2 T 2 ( 2 /)T 2 - 1 (y- 1 )) J R ± T 1± }* 2 

= tr 2 M 2 T 2 {y)T 2 \y- l )T 1± 

= t(y)T 1± . (98) 

Thus, we have proved that the transfer matrix with a particular choice of open bound- 
ary conditions is quantum supergroup SU q (n\m) invariant. 



5 Summary 

In this paper, we have diagonalized the graded vertex model with open boundary condition 
by using the generalized algebraic Bethe ansatz method. In order to get the energy 
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spectrum of 1 -dimensional quantum system defined by equation (23), we assume Vj to 

be zero in equations (58) and (60). However, one would as will assume vf^ ^ 0. In this 
case, equations (56), (69) and (70) lead to the solution of inhomogeneous graded vertex 
model. Formally, one can also define a 1-dimensional quantum system by equation (23). 
Generally, the hamiltonian is not represented in the nearest neighbour interaction form. 
We also show the SU q (m\n) invariance of the quantum spin chain (equivalent to a graded 
vertex model). Thus, the generators of SU q (m\n) commute with the infinite number of 
conserved quantityies. The Hilbert space of the system can be classified according to the 
irreducible representations of SU q {m\n). We hope that it will be help to solve the Bethe 
ansatz equations. 

In order to find the free energy of the system, one should to solve the Bethe ansatz 
equations. Following the method given in reference [18], we can deduce the Bethe ansatz 
equations into those of the periodic case on 2N sites with an additional source factor (see 
ref. [20]). The free energy contains two terms. One is the known bulk free energy, another 
is the surface free energy which is the correction of the open boundary conditions (that 
keeps the quantum group symmetry). This was pointed by de Vega and Gonzalez- Ruiz 
in SU (n) case. 
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6 Appendix A 

The starting point for commutation relations is reflection equation (25). Let a± = a 2 = 
b 2 — 1, b± — b ^ 1, we find: 



A{v)B b {u) = B b {u)A{v) 

ftl2{U+) 11-^21 \ U -)bl 
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i?12(^ + )^2l(^- U6 

' R 12 {u + )\\R 21 {uJ) b b \ 



B b {v)A{u) 



Rl2{u+) ^B c (v)D cb (u) (99) 



Rl2(u + )\i 



Due to eq.(27), it can be checked that the following relation is always true for R(u)\\ 
sin(r] ± u). 



A{v)B b {u) 



Ri 2 (u_)\\R 21 (u + 
R 12 (u + )\\R 21 (u_ 



y/i 

J 61 



B b (u)A(v) 
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R 12 (2u)\\R 12 (u-)\l' 

R ™ {u+) *-B c {v)D cb {u) (100) 



Obviously, commutation relations (29,31) can be obtained from (100). 
Next, let a 2 = 1, a±, bi, b 2 ^ 1, one can get: 

D aibl {u)B b2 {v) = —— — * a B C2 (v)D Cldl (u) 

Rl2{U-) a 1 ll R2l(u+) b 1 ll 

~ R U, Y^ff U, ^i B MD d2b2 {v). (101) 
n 12 {u-) ail K 2 i{u + ) bll 

Substituting (27) to (A. 3), we obtain 

D aibl (u)B b2 (v) = 1 ^ {RM^Z^XSBMD^u) 

Ri 2 {u-) a \ 1 R 21 {u+) b \ l 

- J R 12 ( M _)^i? 21 ( M+ )S 2 1 5 dl ( M ) J D d2fc2 (t;)} 

+ 1 i-jF (102) 

Ri 2 {uJ) a \ 1 R 21 (u + ) b \ l 



fii2(t* + )^iZ 2 i(ti-)gg § g ^^ gc a (t;)A(ti) 



+i2i 2 (u_)S;i22i(«+)Si)*a 1 * 1 ^(«)5 6a (t;) 

+i2i 2 (« + );i;i2 2 i(«_)^ 1 A(t;)S dl («) (103) 



In the following we will calculate the function F for the case of R(u)\\ = a(u) = sin(u+i]). 
The results can be written in a simple form, the main calculation results are as follows. 
Case 1: a\ ^ b\ 

,_, . sinrie l ( u+v * ) sin(2v)sin(u + v)sin(u — v) „ . . ,. . 

^ = <W 2 ( • r , f /• i i -B bl (u)A(v) 

sm{u + v + r])stn{2v + r\) 

sin 2 (ri)sin(u — v) „ . . ~ . . . 

sm(w + v + 1]) 



26 



Case 2: a\ = b\ = b 2 = a 



stn(r])e^ + ^sin(u + v)sin(u — v)sin(2v)sin(r] + e a (2u + r])) 

F = r r r B a [u)A[y) 

sm(2u + r])szn(2v + r])sin(u + v + rj) 

sin 2 (n)sin(n + e a (u-v)) ~ 

—B c [y)U ca \u) 



sin{u + v + rj) 
sin(u — v + r])sin 2 (ri)e^ u ^ 



sin(2u + rj) 

Case 3: d\ = b x ^ b 2 



B c {u)D ca {v) (105) 



sin{rj)e i( - u+v) sin(u + v)sin(u - v)sin(2v)R(2u + r?)? 1 ^ 

F = — — , , —B b Ju)Aiv ) 

sm{2v + r))stn{2u + i])svn{u + v + n) 



sin(u + v + rj) 
sin(u-v + 7])sin 2 (r])e t{u ~ v) ~ 
+ sin(2u + rj) B °M D **W ^ 

Correspondingly, in the case of R(u)™ = sin(r) — u) = w(u),a = l,---,n, the main 
calculation results are presented in the following form. 
Case 1: a\^b\ 

sin(r))e l ( u+v ^sin(2v)sin(v + v)sin(u — v) 

F = b aib2 — — — B bl (u)A(v) 

sinyq — u — v)stn{r] — 2v) 

-K^— 1 ! 1 ^B c (v)D cbl (u) 107 

sin(rj — u — v) 

Case 2: a\ = b\ = b 2 = a 

sin(i])e l( - u+v ^ sin(u + v)sin(2v)sin(2u — 2rj)sin(u — v) 

F = -, : -, : -, : B a (u)A(v) 

sinyq — 2v)szn(rj — 2u)sin{r] — u — v) 

sin 2 {r])sin{r] -u + v) ~ 

B c (v)D ca (u) 

B c (u)D ca (v) (108) 



sin{ji — u — v 
sin 2 (?7)e^"~ 1 ') sin(r] — u + v 



sin{ji — 2u) 
Case 3: a\ = b x ^ b 2 



sin{u + v)sin{u — v)sin{2v)sin(rj)e { ' 4 ) R\ 2 {2u — rj) b l 

F = — — , , r —B b Au)A{v) 

svn{2u — T])sin[r] — 2vjstn{i] — u — v) 
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(109) 



Though we have already simplified the results, they still seem to be too complicated to 
be dealt with. Fortunately, we have found that the results (104-109) can be summarized 
as a concise form which indicates the commutation rules between D aibl (u) and B b2 (v). 
The explicit commutation relations are written in sect. 3. One can prove it by expanding 
relations (30,32) according to different cases mentioned above. Thus, we have obtained the 
commutation relations (29-32). It is also necessary to calculate the commutaion relations 
between B a {u) and B b (v). 

Let ax = a 2 = 1, bi, b 2 ^ 1, we have the results: 



B bl (u)B b2 (v) 



R 12 {u + )\lR 21 {u 



\d1c2 

Kb 2 



R l2 {u„)\\R 2l {u + ) b b \\ 



B C2 (v)B dl (u) 



(110) 



7 Appendix B 

The classical simple graded Lie algebra SU(n\m) is defined by generators h^e^j^i = 
1, • • • , m + n — 1 and the following relations 



[hi, hj] = 0, 

[hi, fj] — a ijfj) [hi, Gj] = —aijCj, 

[fi,ei] = hi,i^m, 

[fmi GmJ-|_ h m , 

[fi,ej] = 0,i^j, 
f 2 = e 2 =0 

Jrn °m w > 
fifi±l ~ 2fifi±lfi + fi±lfi — 

e 2 e i±1 - 2eiei ±1 ei + e i±1 e 2 = (111) 

The last two relations are the so called Serre relations, which are compatible conditions 
for SU(n\m), a^ is the component of the graded Cartan matrix A defined by: 



an = 2,i^m, 

&mm 0, 
a i+l,i = — 1 

a i:i+1 = -l,i^m 

(112) 
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the other elements being equal to zero. We define Oi as: 

Oi = diag(l, 1, • • • , 1, -1, 1, • • • , 1) (113) 

where -1 is th ith element. 

The fundamental representation of the generators takes the form 

id = Eij,i = 1, • • • ,m + n, 

fi = E i!i+1 ,i = l,---,m + n-l, 

€i = E i+hi ,i = 1, • • • ,m + n - 1, 

hi = Wi-w i+ i,i^m, 

h m = w m + w m+1 . (114) 

Here E^ are (m + n) x (m + n) matrices with the element in i-row j-column equal to 1, 
all other elements being zero. 
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